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We introduce a novel model for DNA condensation (whip-toroid transition) using the path integral 
method in the framework of the non- linear sigma model on a line segment. We show that some of 
its classical configurations exhibit toroidal forms, and the system has phase transitions from a whip 
to toroidal phases with a parameter c = S (^r) . We also discuss stability and finite size effect on 
these states. 
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Since in living cells DNA is often packaged tightly, for 
instance, inside phage capsids, DNA condensation has 
drawn much attention [J 0, [E 0, IE IE as well as its 
biochemical/medical importance in the emerging field of 
gene therapy. In fact, when we put condensing agents 
as multivalent cations into DNA solution, it can cause 
DNA to undergo the condensation from a worm-like chain 
(whip) to toroidal states [3, 0, 0, H • 

A double stranded DNA chain can be modelled, for 
example, by a semiflexible homopolymer chain @, IE Q • 
To increase our understanding of "whip-toroid transi- 
tion" , semiflexible homopolymers in a poor solvent con- 
dition (i.e effective interactions between polymer seg- 
ments are attractive) have been investigated as simple 

models M EI El EE EI EE, EE, El 13 EE HE HI 

[221 . I23I |24| . Simulations using Monte Carlo, Langevin 
approaches or Gaussian variational method, calculated 
phase diagram for the semiflexible chain in a poor sol- 
vent 0, [ill, El EE EI EE EE- In theoretical works, 
existing phenomenological models balance the bending 
and surface free energies to estimate toroidal properties 

[HEEEEHEIHIilliEIil. Tt becomes increasingly 
probable that toroid is the stable lowest energy state — 
the ground state. We note, however, that the theoretical 
aspects of the works assume a priori toroidal geometry as 
the stable lowest energy state with no theoretical proof 
[2lj . Moreover, compared to the theory of coil-globule 
transition of flexible chains 0, HE HE HE [2c| w hich 
are well described by field theoretical formalism |2E HE , 
there is no simple microscopic theory, which contains the 
salient physics to demonstrate the whip-toroid transition 
of the semiflexible polymer. Difficulties in formulating 
theory results specifically from the local inextensibility 
constraint of the semiflexible chain 

[El, 

which makes 

the theory non-Gaussian, and also from the non-local 
nature of the attractive interaction along the polymer 
chain, which makes the theory analytically intractable. 

In this letter, we investigate the whip-toroid transition 
of a single semiflexible homopolymer chain with attrac- 
tion, using path integral method and the nonlinear sigma 
model on the line segment. 

In the continuum limit, the Green function (end-to- 
end distribution) of a stiff polymer chain with attractive 



interactions can be given by the path integral 



G{6,R;ui,u f ;L,W) = J\f~ 



r(L)=R,u(L)=u f 

V[r(s)} , 

r(0)=0,u(0)=3; 



-H{r,u,W] 



(1) 



with a constraint |u| 2 = 1 [E HI- s is the proper time 
along the stiff polymer chain of total contour length L. 
f(s) denotes the pointing vector at the 'time' s in our 
three dimensional space while u(s) = corresponds 
to the unit orientation vector at s. Af is the normali- 
sation constant. Following Freed et al. (EHE|, the di- 
mensionless Hamiltonian can be written by TL[f, u, W] = 

J ds [H(s) + Vat(s)] where H(s) and Vat(s) are the lo- 
cal free Hamiltonian and the attractive interaction term, 
respectively: 



H(s) = - 
V ; 2 



d 
as 



(2) 



V A t(s) = -W I ds'S(r(s)-r(s')). (3) 
Jo 

I is the persistence length and W is a positive coupling 
constant of the attractive interaction between polymer 
segments. Thermodynamic (J — X/iksT) is implicitly 
included in I and W . This will be revived when we con- 
sider the thermodynamic behaviours of the system. In 
what follows, we express r by the unit tangent vector u 
and therefore the Hamiltonian Ti(u) in terms of u. Hence, 

the Green function G [vl, R; Ui, Uf \ L, W^j becomes a path 

integral over u with a positive coupling constant W , re- 
gardless of r, 



G 



V[u{s)\ 5 ( Jq ds u(s) - htj e 



-H[u,W] 



(4) 



where we used r{L) = j ds u(s) and the Jacobian is 
absorbed by M which is neglected here. The delta func- 
tion selects out the end-to-end vector. When W = 0, 
our free dimensionless Hamiltonian with the constraint 
l^( s )| 2 = 1 can be interpreted as the low energy limit 
of a linear sigma model on a line segment, or quantum 
equivalently a nonlinear sigma model on a line segment, 
rather than some constrained Hamiltonian system. 
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We begin with (9(3) nonlinear sigma model on a line 
segment which is nothing but a quantum mechanics of a 
limited time s € [0, L\ with a constraint. The constraint 
|u| 2 = 1 restricts the value of u on a unit sphere S 2 . Its 
polar coordinate decomposition in the free Hamiltonian 
7i(u) = | J Q L ds \du(s)\ 2 gives the action: 

S[6»,<p u ] = l -J ds [(de u ) 2 + sin 2 e u (d^ u ) 2 ] . (5) 

This is called the nonlinear sigma model since the action 
is 0(3) symmetric but some of its transformations are 
realised nonlinearly. Minimizing the action {5| in terms 
of 9 U and ip u yields the classical equations of motion: 



sin2# u ,„ . o 
-d 2 + ^^(dpu) 2 



d 2 + 2(d9 u )cot6 u d] <p u 



(6) 



Our aim is now to explore classical solutions of eq.© 
and to study the whip-toroid phase transition in the pres- 
ence of attractive interactions. Consider classical solu- 
tions of eq.© with a trial solution 9 U = 0. The first 
equation of © leads to sin 29 u (ip u ) 2 = 0. Thus the so- 



lution is either 9,, 



0, 77, 7r or tp v 



0. The solutions 



9 U = 0, 7r or ip u — with 9 U = are equivalent to having 
a constant u. Accordingly, classical solutions reduce to 
9 U = | or m = const. Substitution of 8 U = \ into the 
second equation of motion (J6j> gives d 2 tp u = 0. Therefore, 
we have two types of classical solutions: 



u(s) = const. 



= — and ip u = as 



(7) 



where a,b are constants. By symmetry argument, we 
state that the solutions (J7J) represent the general solutions 
[29l | . That is either a constant u(s) (rod solution) or a 
rotation at a constant speed along a great circle (toroid 
solution) . 

Now we consider the attractive interaction term ((3|). 
It is difficult to interpret it in the context of quan- 
tum theory due to its non-local nature along the poly- 
mer chain. However, we can solve them with our clas- 
sical solutions (J7]). Let us rewrite eq.© by r(s) — 

r(s') = Jq dtu(t) — Jq dtu(t) = J 3 ,dtu(t), that is, 
Vat{s) — —W Jq ds' 6 (J* dtu(t)). Hence the problem 
is now reduced to the one in the u space: finding non- 
zero values of 5 (J s s , dtu(t)^ with the classical solutions 
([7]). That is to find u(s') for a given s, which satisfies 
J^, dtu = . In the polar coordinates, this is expressed 
by J^, dt sin 9 U cos (p u — 0, J^, dt s'm9 u sin ip u = 0, and 
J s , dt cos 9 U = 0. The first classical solution does not 
satisfy these equations and thus derives no attractive 
interactions. If we substitute the second classical solu- 
tion of eq.JT]) into the equations, we have cos9 u (s) = 0, 
J 3 , dt cos(at + b) = i (sin(as + b) — sin(as' + bj) = and 
J 3 , dt sin(at + b) = i (cos(as' + b) — cos(as + b)) = 0. 
Hence we have solutions: s — s' — 2nn/a > 0, n G Z. 



Without any loss of generality, we assume a > and 
n e Z_|_. Introducing N(s) = [as/2Tr] by Gauss' symbol 
P, we obtain J^_ 2v/a dtu(t) = f*_ iw/a dtu{t) = • • • = 
fs-2TrN(s)/a,dtu(t) = 0. Therefore, the attractive poten- 
tial is given by V AT (s) — —W ■ N(s). Note that N(L) 
represents the winding number of the classical solution 
([7]) along a great circle of S 2 . Finally, an integration over 
s yields the dimensionless Hamiltonian with our classical 
solutions: 



U [u,W}= [ dsH(s)+ [ dsV AT {s) 
Jo Jo 



LI , 
- ~2 ~ 

LI , 
= — a 2 -WL 



2k sr^ , 2ir I aL „ T/rN \„ T/rN 

— E k+ — U- _iV i P L 

n ^ — ^ n \ /.it I 



k=l 



N(L){l-^(N(L) + l)} 



(8) 



The first term denotes the bending energy, and the sec- 
ond and the third terms are thought of as 'topological' 
terms from the winding number. 

The non-zero winding number of the classical solution 
in the u space means that the polymer chain winds in 
the r space as well. That is, when a >j-, configurations 
around the second classical solution ([jj start forming a 
toroidal shape since 



i {sin(as + b) — sin(&)} 
r(s)— I — i {cos(as + b) — cos(6)} 
const. 



(9) 



and stabilise itself by attracting neighbouring segments. 
We call such classical solutions the "toroid states." 
Whenever a increases and passes through the point ^il 
for n £ Z_|_, another toroid state appears with the in- 
creased winding number n. Note that the radius of the 
toroid state is given by -. When the chain of contour 
length L winds N(L) times we have the N(L) circles of 
each length ^ and the rest (L - ^-N(L)). The second 
and third terms in the second line of eq.© result from 
the former and the latter respectively. When < a < =j-, 
the chain cannot wind like the toroid states. Both ends of 
the chain are not connected to each other, thus can move 
freely as well as any other parts of the chain fluctuate. As 
long as the total energy of the chain does not exceed the 



bending energy of 



at a 



2jl 

L 



, they can whip with 



zero winding number. We call such low-energy states 
the "whip states." Although the definition includes fluc- 
tuations around the classical solutions, unless otherwise 
stated, we primarily refer to the classical solutions of such 
states, which are rather bowstrings than whips. 

We now explore the exact energy levels of the whip 
and toroid states, and discuss the phase transitions be- 
tween these states. The dimensionless Hamiltonian of 
the second classical solution (J7J) is a function of l,L,W 
and a: 

H cl (a,l,L,W) 

H a2 + irW N ^ ( N ( L -) + 1 y WL . N ( L ^ ( 10 ) 
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This matches with the first classical solution when 
N ^ = for a = is defined. Accordingly, the above 
expression is valid for all classical solutions. 

Consider first a case with L, W, and I fixed. By defi- 
nition, Ti(a) = Ti c i(a, I, L, W) is continuous in the entire 
region of a > and is a smooth function in each segment: 

~27rjV 2tt(N+1) ~ 
L ' L 



a G 



for JV £ Z; 



However, it is not 



smooth at each joint of the segments: ^ G Z+. We plot, 
in FigQ] the energy level as a function of a for different 
values of c, showing qualitative agreement with Conwell 
et al. ||. Given N(L) = N is fixed, the Hamiltonian fit)]) 
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FIG. 1: The dependence of the energy Tt(a) on x — aL/2n 
and c. Ti.(a) is scaled by the factor of yWl for convenience. 
The definition of c will appear shortly. 



takes a minimum at a = a c (N) = {^-N{N + 1)) 1/3 . 
Accordingly, each segment falls into one of three cases: 

(i) When a c (N) < ^7—, H(a) is a monotonic function 
in the segment and takes its minimum at a — ^j^- . 

(ii) when < a c (N) < 2 " ( ^ +1) , H(a) behaves 
quadratic in a and takes its minimum at a — a c (N). 

(iii) When 27r ( j ^+ 1 ) < a c (N), TC(a) is monotonic in the 



segment and takes its minimum at a 



27r(Af+l) 



The first and third cases are physically less relevant since 
they mean no (meta-) stable point in the segment. So, we 
focus on the second case. Introducing a new parameter 

c = (j^) the condition for N or a to be in the second 
case turns out to be 

N L {c) < N <Nu(c) forc>4, 

I < N < N v {c) for < c < 4, (11) 



the whip phase to whip-toroid co-existence phase transi- 
tion may occur. On the other hand, when c > |, there 
always exists at least one (meta-) stable state for some 
non-zero value of a with positive winding number N(L). 
The number of minima are roughly given by the width 
of the region for N, i.e., Nu(c) — Nl(c). For example, 
when c > 4, Nu(c) — Nl(c) > 3, and therefore there 
are at least three minima with positive winding numbers 
greater than 1. When < c < 4, the condition of hav- 
ing three minima is c > f . To summarise, when c > § 
there exist at least three minima with different winding 
numbers. This occurence of the multiple local minima is 
because the energy is given by the balance between the 
bending energy and the attractive potential energy. The 
former is a monotonically increasing function of a while 
the latter is monotonically decreasing but not smooth 
function of a. This non-smoothness is a source of mul- 
tiple stable states. Note that our precise analysis shows 
that the number of minima could be reduced in some 
cases, for example, one can confirm it shortly by plot- 
ting the Hamiltonian (fT!?]) with some finite size effect on. 
The existence of multiple minima indicates the first or- 
der phase transitions between these stable states. In fact, 
when the energy of the N = 1 stable toroid state becomes 
zero, the whip dominant to toroid dominant phase transi- 
tion occurs. Such a value of c is 27/16. When c > 27/16, 
the toroid states will dominate the action. Further dis- 
cussions on the phase transitions can be seen by changing 
the value of c [30]. 

One can plot the critical value of c where the minimum 
of the N-th. segment emerges and vanishes. The lower 



bound of the N-th segment is c 



(N)_ 



bound is c 



(iV)_ (N+iy 



N+ 



j , while the upper 



. So, when < c < c[f\ the A^- 



th segment has a minimal and (meta-) stable point. For 
example, when i < c < 4, the first segment a G [^r-, 
(i.e. N — 1) has a minimal point at a — a c (l). 

So far we have dealt with the classical solutions, which 
are derived from the first derivative of the action. Thus, 
they may correspond to the global/local minima of the 
action in the configuration space. However, the solutions 
are not necessarily stable unless we take into account the 
attraction, since the second derivative of the action is 
zero or even negative in many cases. When we take the 
attraction into account, we observe the followings. When 
c > 4, the toroid states become stable under the quantum 
fluctuations, since deviations from such toroid states will 
cost a large penalty in energy. When c is small, a number 
of the whip states become equally or more probable to 
any toroid state. 

For large c, the ground state — the dominant toroid 
state of the winding number A c can be estimated by the 



, , . inequality relation of c: < c < c\f^ . It reads A c ~ c. 

where N = [aL/2n], N L (c) = §(l-§ + \/l _ c)' and Using this, we can estimate that the radius of our ideal 



Nu(c) = % (l + yiTi). 

When Nu(c) < 1 (i.e. c < 5), the above second con- 
dition vanishes and thus the whip states only survive at 
low energy. Therefore, at the critical value of c = |, 



toroid behaves 



2ttAT c 



4ttI 
WL ■ 



Our ideal toroids have 



zero thickness, but the finite size effect of their cross sec- 
tions can be approximated by the hexagonally arranged 
DNA chains with a van der Waals type interaction, i.e., 
with the effective nearest neighbour interactions. In the 
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case of N(L) > 4, it leads to the modified Hamiltonian: 

H(a) = ^ - 2 -^VD(N(L)) 
2ttW { aL 



a \ 2ir 



N(L))Gap(N(L)), (12) 



where VD(N) = 3N - 2y/3yfiST - 1/4 and Gap(N) = 
VD(N + 1 ) — VD(N). Following the same procedure for 

< c < , one can find N c ~ (2\/3c) 5 for large 
c. By r c 



— L 



. we now estimate the mean radius of 
the toroid (i.e. the average of inner and outer radii) in a 
physical system in more detail. A coupling constant can 

be given by W = where k is the number of the 

electric dipoles in a monomer segment, which create van 
der Waals interactions of the magnitude e. l m denotes the 
length of the monomer along the chain contour, taken to 
be l m ~ 5 bp — 1.66 nm in the end. Substituting N c and 
the above, we obtain 



r, - (6^L* ( i-Y = (6tt)-*£*(U0 S k< 



W 



knT 



This L= dependence agrees with the proposed exponent 
in the asymptotic limit pj, S2, [H, [2| . 

We estimate the mean toroidal radius of T4 DNA in 
low ionic conditions reported in ref.[4j. Using L — 57/im, 
Z ~ 50 ~ 60 nm, and l m , the mean radius of the toroid is 
r c = 29.09B-§ - 31.29S-S [nm], where B = -^ T . This 
result is in good agreement with the experiment r c ~ 



28.5 nm for i? ^ 1.15. The same argument for Sperm 
DNA packaged by protamines (L — 20.4/xm) gives the 
analytic value r c — 23.695" ~ 25.48-B" [nm], which is 
comparable to an experimental result r c ~26.25 nm. Note 
that the latter has a larger diameter of the segment and 
is expected to have the weaker interaction with smaller 
B. 

The exponents v predicted in the literature for r c ~ L v 
are v — i in most cases [3 [22|, HU, [24[. However, 
they are inconsistent with the experimentally well known 
observation that the radius is independent of the chain 
length [2|, y, |j] . This might suggest that the real interac- 
tion is not van der Waals like, or at least is not a single 
van der Waals type interaction. It should be noted here 
that combinations of our ideal toroid and its finite size 



An- 



effect can give a range of v = — 1 ~ ^ in some region, 
other important and interesting remark is that, in fact, 
when we apply Coulomb like interactions to our approx- 
imation, we observe some asymptotic behaviour that the 
radius remains constant as L changes. The precise anal- 
ysis is to be studied and given in the near future. 
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